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AN INVESTIGATION OF ASYMPTOTIC 
SOLUTIONS OF ORDINARY DIFFERENTIAL 
EQUATIONS 

By Robert M Myers 


SUMMARY 

This paper presents a brief study of power series, Frobenius 
and Thome 1 solutions of ordinary differential equations In 
addition, an elementary proof of the existence of actual 
solutions corresponding to the formal Thome 1 solutions is 
given Moreover, a numerical method is developed which 
enables one to continue the accurate values of the Thome 1 
solutions to small values of the independent variable 

INTRODUCTION 

Consider the following N th order homogenous system of 
ordinary differential equations (ODE) 


Cl) 


dy (x) 
~dX 


xh lX x 

n=o 


-in 


y(x) 


where 


CO 



n=o 


converges for 


integer 


xj > a > 0 and h is an 


1 



It is quite natural to attempt to determine solutions of 
the form 


CO 

(2) zoo - XX x ' n 

n=o 

If (2) is formally substituted into the system of 0 D E ’ s 
(1) one will obtain the equation 


(3) 


CO 



n=2 


(n-l)x n _, x 


- n 



X 


-n 


Assume for now that A is a nonsingular matrix and 
h > 0 Then, one obtains by equating coefficients of 
like powers of x 

ll — 

x £ " A 0 X 0 > Zq = £ since A q is nonsingular 

(4) x h_1 . 0 = A y + A y = A y. ->■ y_ ~ 0 

— l^-o o*~l o—l ^-1 — 


1 = A 2X, + Vl + A o— 2 = A oi 2 - 0 


• 

• 

• 

« • • 

• ♦ 

x° 

0 = 

A hXo + 

• • • + A v + 

0 

i 4 

ii 

o 

< 

X - 1 

0 = 

Vllo 

+ + A v + 

lZh 

A oZ h+1 - Z h+ 1 

x- 2 

= 

A h*2t> 

* ’ ‘ ' + A A+i 

+ V-h + 2 

but 

*1 = 

o - Z h 

« - i 
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It xs clear that substitution of a power series solution 
into the ODE when A q is nonsingular and h > o leads 
only to the trivial solution, 1 e , one obtains only the 
solution y(x) = J3 Also it is rather exceptional to 
obtain nontrivial power series solutions for the case 
where A q is singular 


Now consider the case h = -1 and A does not have a 

o 

nonpositive integer eigenvalue Once again, one can attempt 
to substitute formally a powei series solution into the 
differential equation 


( 5 ) 


xdy(x) 

dx 



yW 


CO 





-n 


One can attempt to solve for y by equaling like powers 
of X 


x° 0 = A y y =6 since zero is not an 

— o—o *—o — 

eigenvalue of A j 

(6) ° l 

x -1 -y, = A„y + A y, = Ay, but -1 1 is not an 

^-1 1^*0 o—l o-^-l i 

eigenvalue of A \ y, = 0 

0 o *-l 


Similarly, one can conclude X 2 ’— 3’*** are F^ 1 ' 0 vectors 

Hence, the formal procedure leads only to the trivial 
solution y(x) = 0_ 


3 



Now consider the case h 


-2 


(7Ha) 


CO 


Xr (“'‘'In-l 1 

n-2 


00 ! n 

^E^ A n-2-kXk 

n=2 \k=o 


n 


(b) 



choose an arbitrary value of y_ o and 
this determines y 



-2y = Aj + A ^ 

"2 l^o 0^-1 

clearly y 2 is determined from 

y and y, 

— o —1 


- (n+1 ) 

i (\ 


-ny = A _y A .y n + * • • 

^-n n-1— o n-2^1 

+ A y , thus, y is determined by 
o— n-l — n 


A 

n- 


1 






1 


Notice that the determination of y^ proceeds regardless of 
the exceptional properties of the lead coefficient matrix 
A^ Indeed, the power series solution can be used when 
h < -2 and h is an integer So far, the procedure out- 
lined above has been strictly formal, l e , no justification 
of the validity of the formal solution has been given What 
is required is a justification of the procedure of term-by- 
term differentiation and equating like powers of x It 
is somewhat surprising that this procedure is completely 
justified provided one establishes convergence of the formal 
pouer series (see [1] or [3]) For proof of conveigcnce, 
see [2] or [4] 
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Now, return to the case h = -1 For simplicity, consider 
N = 1 , i e , the order of the system is one (Scalai problem) 


( 8 ) 


xy * 



y 


For this problem, a slightly more general substitution turns 
out to be successful Consider 


( 9 ) 



CO 



n=o 


Substitute (9) into (8) and equate like powers of x 


( 10 ) 



(Clearly this equation will be 

satisfied for a = a and 

o 

arbitrary y ) - 

x a_1 (a o - l) yi = a 0 ^x * Vo 

or -?i = a Vo 

(Thus , once a value of y Q is 
assigned, then y^ is determined) 


< a o • 

- n 

)y = a y 

+ Vn-l 

+ • 1 

• + 

a y -ny 

= Vn-l 

+ • * ' 

► -f 

a y 




n o 



IS 

determined 

from 

a , 

n 


’ * * ’ ^n- 

P 


5 



The solution (9) is referred to as a Frobemus solution 
The more general problem of an Nth order system also admits 
solutions of this form There are certain exceptional cases 
that must be dealt with separately for N > 2 Since there 
are several good references on the so-called regular 
singularity [2] or singularity of the first kind [4], the 
case h = -3 will not be treated more fully here The 
formal manipulations can be a ■postev'Lov'i justified by merely 
establishing convergence of the expansion (see [4]) 

Now, consider the case h > 0 and N = 1 


(ID 


y'(x) - fx h y^a n x n )y 


n=o 


Attempt to determine a solution of the foim 


( 12 ) 


y(x) = I >^y n x n ! xC 'e i?(x) 

,n=o 


where 


h+l 


P( x ) = p h+l X + *’* + p i X 


6 



If (12) is formally substituted into (11), then one obtains 


+ X> n * a - n (o' + np h+1 -' 

(n=o n=o 


(13) ♦ •• ♦ 2p 2 x ♦ Pl j|e*«*> - I *J2 J = n -' ,, U 

,n=o 


a-n p(x) 

y x 
n 


,n=o 


Formally solving (13) leads to the relations 


h+ 1 h+1 


(14) 


Pi = a 


h 


a = a 


h+1 


y Q is arbitrary 


' y l ~ a h+2 y o 


- 2y 2 = a h+3 y o + a h+2 y l 


_ny n = a h+n+i y o " a h+n y i 
+ -* +a h+2 y n-l 
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The success of formal solutions (Thome 1 solutions) of the 
form (12) have been known for many years Unfortunately, 
these solutions seldom lead to convergent expansions when 
N > 2 Thus, for this case (h > 0) one can rarely justify 
a yostev'iopv the formal manipulations Around the turn of 
the century, Poincare' proved that there is a connection 
between actual solutions of the differential equation and 
these formal Thome' solutions 

The integer h + 1 is referred to as the rank of the 
singularity at x = « 

THEORETICAL INVESTIGATION OF THOME' SOLUTIONS 

It is best now to consider a particular problem If the 
problem is carefully selected, the general theory will be 
illustrated but perhaps the involved analysis can be 
simplified In any event, one would wish to develop an 
intuitive understanding of the general theory and perhaps 
eventually read one of the excellent references ([4], [6] 
or [8]) available to more fully appreciate and utilize the 
general theory 

Consider the third-order system with a singularity of rank 
one 


(15) (x) - | diagonal (0 , 1 , -1) + i diagonal ^-1 ,^ 7 ) 


+ -j diagonal (2 ,l,10)jy(x) 


S 



(16) let 


A 

o 



A 


2 


A 

n 


00 


diagonal (0 , 1 , - 1) 


diagonal -1 


II) 

’ 2 ’ 2 / 


diagonal (2,1,10) 


0 for n 2 3 


A + — A + A. 
o X 1 2 2 

X 


Since the coefficient matrix A(x) is diagonal, one can 
immediately determine a fundamental solution (nonsingular 
matrix solution) See [4] 


(17) Y(x) 

51,00 

y 2 (x) 

Z 3 Cx) 


= z, (XJ X-2 OO >Z 3 (x) 



1J 


Now, this problem is much too simple to illustrate any 
general principles Theiefore, the problem is transformed 
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(.18) let 

where 


v(x) = POOz(x) 
P (x) = I + ~ C 


I = diagonal (1 , 1 ,1) 


C = 


1 

1 

3 


2 3 

A 9 
1 1 . 


It is easily sho\% r n 


(19) 

£' (x) 

(20) let 

B(x) 

then 

PB 


that y(x) satisfies (19) 

= {P -1 (x) (A(x)P (x) - P'(x))}£(x) 

= P -1 (x) (A(x)P(x) - P'(x)) 

= AP - P' 


(21) assume B(x) 


00 



The relationships (22) are a result of substituting (21) 
into (20) 


( 22 ) 


B 

o 



for n > 5 



B 

n 


A 

o 


A + A C - CA 
2 o o 

A C - CA 2C 


a 2 c 


CB 


2 


-CB 


n- 


2 
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.*.(23) 

v ’ (x) 

- 

{/ 

} ^ 
V 


cc 



Since 

yy x- n 

n 

is 

small 


r_=2 



to compare solutions 

£0) 

(24) 

w * 

(X) 

_ ( 
“ i 

Recall 

that B 

O 

A 

o 

and 


B, 


x \ \ 


x / 


x.A x " u (x) 

n-2 ) 


h = A i ’ 1 e * -o -1 


B and B. 
o 

axe diagonal matrices and hence one can immediately obtain 
a fundamental solution W for (24) 


(25) 


W(x) = [w 1 (x),w 2 (x),w 3 (x)j 


-i (x) = 


1] x 

oj A 


^ 2 (x) = 


~i 1 


2 x 
x e 


yt 3 00 


rni 1 


2 -x 
x e 


One might ask what connection, if any, exists between tiie 
solutions of (23) and (2^) There is a remarkable connec- 
tion first established by Poincare* and -ex tended by many 
other researchers (see [8]) 
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THEOREM If j arg x| < > there exist actual solutions 

ZlCx) .£ 2 ^ x ) ’I 3 of C 23 ) 

-» 

such that 1, MXiOO - w ± (x)!i = llvjl - °( 

i I £,00 - w 2 fx)| | = I lnJ I ^ o( 

ill; ! I £3 00 - w 3 (xj I I = i I w 3 M x o( 

where j j Z | | = | j 

implies (in the context of this 
problem) that there exists a constant M , such 
that, 0 < M < 00 and. |f(x)j < — as x -> 00 » 

Notice that [ |£ 2 (x) " 1-2 00 1 i 1S not required to tend to 
zero Since ||w 2 0 O|! is exponentially increasing the 

need not tend to zero However , 

the relative difference becomes small, 1 e , the number of 
significant figure agreement betveen £ 2 (x) and w_ 2 (x) 
increases with increasing x The reader may be familiar 
with Stirling's formula for n 1 Precisely, the same 
behaviour is observed, 1 e , Stirling's formula yields more 
significant figure agreement with n* as n increases 
The difference between n' and the Stirling approximation 
is unbounded as n rends- to infinity 


product, | jw 2 (x) | | x 0^; 


and f (x) 


0, l 


'll 

^2i 

J, 


-II = max lyj 
1 < 1 < 3 1 
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r-H I X r- i I X ' — I l X 



Now, return to the problem of establishing the existence of 
actual solutions y^(x) 5 such that [[^(x) - w^(x)|j 

! |w x oo 1 1 x 0^ 

co 

(2(5) let E (x) = 2^B n x" n 

n =2 

then (23) becomes 

(Z7) £'(x) = {b 0 + | B L + E(x^JxCx) 

One can also introduce a norm for a matrix A possible 
compatible norm is given m (28) 


(28) | [ E | 


I I [e 1=! ] | [ 


max 

1 < 1 < 3 



It is easily shown that for this choice of the norm of a 
matrix, the usual norm properties (29) are satisfied 


I IB*] | * l|E|| x ||*|| 

ilE^M < IJEJI x | |E 2 | I 
Me, + e 2 I| < | |E 1 1 ] + ||b 2 || 
j j E | j = 0 IFF E = 0 

i j aE M = | a j x j 1 E | J 
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one might 


Since E (x) (2 6) 

expect E(x)y_(x) 

( B o + I B i)lM 

( B « + h B 1 )i« > 


is small ^ | | E (x) | | = 0 

to be small m comparison 
In fact m comparison wi 



with 

th 


th® term E(x)y(x) is approximately 


known, 1 e , this term is approximately zero Thus, (27) 
takes the form of an mhomogenous equation where E(x)y_(x) 
is treated like an mhomogenous term 


One can attempt to use the method of variation of parameters 
to obtain the general solution of (27) If E(x)£(x) is 
precisely known, then one can obtain the general solution 
of (27) m the form of an indefinite integral Since £(x) 
is not known, one arrives at an indefinite integral equation 
for the general solution of (27) 

The homogenous portion of (27) is merely (24) Recall that 
a fundamental solution of (24) is given by (25) 

Assume a solution y.(x) of the form 

(30) £(x) = W (x)c (x) 

where £(x) is to be determined Substitute ( 3,0 ) into (27) 
then one obtains 

(31) y(x) = W(x) E(s)y(s) ds 


14 



It is easily shown that any continuous solution of the 
integral equation (31) is a differentiable solution of (27) 
This is one big bonus that integral equations possess and 
it can be expected that the analysis can proceed more 
easily if one considers the integral equation (31) rather 
than the differential equation (27) 

OBSERVATIONS (See (25)) 

1 w L (x) dominates w (x) for x increasing to infinity 

w (x) dominates w_ 2 (x) for x starting at large values 

and decreasing to small values of x . (Recall 

|arg x| < y) 

2 w 3 (x) dominates w^ (x) and w 2 (x) for x decreasing 

from a large x value to a small x value 

3 w_ 2 (x) dominates w^ (x) and w (x) for x increasing 

to infinity 

Intuitively, what is meant by the dominance of one solution 
over anothei is the relative growth of one solution compared 
to another For example, w^Cx) dominates w 2 (x) for x 

decreasing For decreasing x , w (x) is an increasing 

function and w (x) decreases exponentia] ly fast Hence, 
w x (x) dominates w 2 ( x ) for decreasing x 

One can use the above observations to construct integral 
equations which force the dominance of a solution which 
behaves like w (x) Three different integral equations 
will result, each corresponding to the three distinct 
solutions w^(x)(i = 1,2,3) 
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Recall W(x) 


diagonal 



(32) 

(a) 


W _ 1 (s) - 

diagonal ys ,s 2 


(b) 

let 

W 1 (s) - 

diagonal^O , 0 , s 


(c) 

let 

w 2 cn - 

( ‘I 

diagonal! s,s 

clearly 


(s) + 

lf 2 fs) - 

W“ 1 (s) 

(33) 

(a) 

let 


= W(x)W 1 (s)E(s) 


(b) 

let 

K 2 (x ,s) 

= W(x)W 2 (s)E(sl 


-s 


1 

2 


“S 



One should note that W(x)W (s) includes the effects of 
only W 3 (x) which is the solution of (24) which is dominated 
by w x (x) for increasing x W(x)W (s) includes the 
effects of vn^x) and w 2 (x) Recall that w_ 2 (x) 1S 
dominated by w^Cx) for decreasing values of x x 


Consider the following integral equation 


(34) l x (x) - ^(x) + K 1 (x,s)y 1 (s) 


ds 


+ K 2 (x J s)y 1 (s) 


ds 
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For simplicity , restrict the remaining portion of the 
discussion to the case arg x - 0 , 1 e , x is real and 

positive 

If a <. x < «> , then the integral from a to x m (34) 
proceeds with ds > 0 and the integral from =° to x 
proceeds with ds < 0 Due to the above split that has 

been performed (see (32) and (33)), one has forced the 
dominance of a solution of (27) which behaves like w^ (x) 
which is a solution of (24) It is relatively simple to 

t 

differentiate formally the integral equation (34) to 
demonstrate that continuous, absolutely mtegrable solutions 
of the integral equation are indeed differentiable solutions 
of (27) 


The proof of the theorem now proceeds m two stages First, 
it is necessary to establish the existence of a solution 
2~(x) to (34) Secondly, it is demonstrated that 

Ihpx) - w x Cx)|| = I |w x (x) | | X o(|) 

In order to establish the existence of a solution to the 
integral equation (34) , consider the Picard iteration 



(3S) 


£“ +1 (x) = w x (x) + f K 1 (x,s)2”(s) 


ds 


+ | K 2 (x,s)y n (s) ds 
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Since 

1 |E(S)|| 

= o(y let 

M be chosen so that 

(36) 


1 |B(s) I I 

< 

- 2 
s 


LEMMA 

I£ & > max(4,6M) and 

a <' x < oo 

then 

(37) 

iiz^w ‘ 

£>)ll i (§) 

1 |w x (x) I I 

■ (i) n - 

PROOF 


HM - 

_0 




ij; (x) = 

5L-,. (x) 



1 lijcx) 

- £°C*)II = 

1 1—3.1 1 ‘ 

(i)° Hi," 


*. fox n 
some n 

= 0 the lemma is true 
- 1 and attempt to show 

Assume true 
true for n 


II£" +1 M - 2 ;cx)|| (|)V (X - S » jL llgcj 


-ss^wii i d *i * (i)^ii2.;w 

*co S 

- £i _ 1 C s ) II l ds I 
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X|)-> 



Fiom the inductive hypothesis one concludes 





- (x-s) 


< 1 fox 


a > 4 



3M 


x 


2 


6M < a < x 
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For a > 4 


it is easily shown that 



monotonically increasing with s for fixed x on the 


20 



interval 


a <. s < x 


It is also easily shown that 



for 4 <_ a <_ x 


Thus , 









M 

3 

s 



ds 


< 


✓2 

a 




(|) n x I ISxCx)! I 




M^QO I I 


The results of case 1 can "be applied to obtain estimates 


for 



i 


i 



4 - 




X 



or 


llxfhx) - II 
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It is readily verified that 


( 38 ) 


~n+l 

*1 


00 


n 



lc=o 


Therefore, the above lemma establishes the uniform conver- 
gence of the above sum as n -*• 00 and hence (x) 

converges uniformly Denote this limit vector by y (x) 

2^00 is clearly a continuous solution of the integral 
equation In addition, the above telescoping sum coupled 
■with the triangle inequality yields an estimate of the 
growth of | |y*(x) | | 

(a) I I y 1 + 1 Cx 3 I I < 2(^) k I I | 

k=o 

( 39 ) 

o>) \ ll^ooll <2 1 1 c^) | j = | 

This estimate of the rate of growth of y^fx) 1S sufficient 

to establish that 11x^00 ~ 00 | I < ||w 1 (x)|[ x 0 
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(40) | ly^Cx) - w L (x) | 


ds x n e' 


3M , 2M 


_ x / 2 

Clearly e is an exponentially decreasing function and 

hence, the first term on the right-hand side of (40) is 



The entire expression is 



or 


| | w x (x) 1 | x 0 This establishes the existence of an 

actual solution of (27) which behaves like w (x) as 
x ->• 00 . 


One could construct other integral equations to prove that 
there exists an actual solution y (x) of (27) which behaves 
like w_ 2 (x) as x -*■ 00 , and similarly one could demonstrate 
a connection between a y (x) solution of (27) and w_ 3 (x) 
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This method of proof is not going to be repeated for y (x) 
and y 3 (x) Only the integral equations will be determined 


Recall that 
increasing 

Z 2 (x) 

X 

dominates w_ x (x) and w (x) 

Thus, let 

W 3 (x) 

/ - X/ 2 s\ 

= diagonal i s , 0 , s e ) 


W 2 (s) 

= diagonal ^0,s 2 e s ,oj 


K 1 (x,s) 

= W(x)W 1 (s)E(s) 


K 2 (x,s) 

= W(x)W 2 (s)E (s) 


( 41 ) 


h (x) 


w 2 (x) + 



K 


l(x ; 


0z 2 cs) 


ds 


+ 


x 


K ,( x >s)l,(s) 


ds 


Since w_ 3 (x) do6s not dominate any solution for increasing 
x , the corresponding y_ 3 (x) integral equation is 


(42) 


Z, (x) 


w,(x) + 


W(x)W 


-I 


(s)E(s)y 3 (s) ds 


So far it has been possible to avoid using Thome’ solutions 
m an attempt to establish a connection between (27) and 
(24) 
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It turns out that for the particular problem under investi- 
gation, one is able to determine (guess) by inspection the 
first few terms m the Thome' expansion 

Now, consider a different method of attack which more closely 
parallels the approach developed earlier m this paper That 
is, attempt to determine solutions via formal substitutions 
Since, the rank of the singularity of (23) or (27) is one, 
the formal solution takes the form 


(43) 


y(x) 



a 

x e 


Xx 


If one formally substitutes (43) into (23), then one obtains 

(44) Ay = B y 

*-o o^-o 


One clearly obtains the values 0, 1 and -1 for X and 

for the corresponding- y^ vectors 
Actually any multiple of these vectors is satisfactory 


"1" 


“0" 


"O’ 

0 

0 

> 

1 

0 

and 

0 

1 


Corresponding to any of the three possibilities for X , 
one can formally determine a ,y 2 ,y_^ , • • • For example, 

consider X = 0 The corresponding formal solution is 
of the form 


(45) 


CO 



11 = 0 
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Vi + B lEo 


a = -1 

(46) -2 £l = B q £ 2 + B lEl + B 2Eq 

-3e 2 - b oE3 + b iE2 + b 2Ei + b 3Eo 

' 4 H 3 = V4 + B l &3 + V2 + Vl + Vo' 


If one solves the relations (46), then one obtains for 



Note that w^x) = P. 0 x \xj w bich is merely the result 

of truncating the formal sum (45) after the first term 


One could proceed to recursively determine 2 4 >E 5 j*** • 

It turns out that if one considers the truncated formal 
% ' 

E (Y - V 

p^x , then there exists an actual solution 

k=o 

£ 

y 1 (x) of (24) such that My^U) -^P k x or “ n || = ** X 0 

k=o 

This result was established m the previous theorem for 
i = 0 This more general result will not be established 
here (see [4] , [8] ) 



26 



Poincare’ introduced the folloi^ing definition A function 
f (x) is said to be asymptotic to 

N 


a x n as x m 

n 


n=o 


the sector / if |f(x) 


property is denoted by 


v' k i ’ °(:i4r 


This 


k=o 


(48) 


00 

'v ^ 


fCX) ~ 2_A X '" > X 6 ! 

n=o 


3 

If one evaluates y~^p n x~ (n+1) at x = 100 one obtains 

n=o 


(49) 


' 009801007’ 

- 00000098 

- 0000029 4_ 


Since one can obtain exact solutions via the transformation, 
(18) one can easily s how that 


(50) 


3 

| INEXACT - ][> n x- n | 
n=o 


3 x 10~ 9 


Since | [y (100) || & 01 , one thus concludes that the formal 

1 

truncated solution y 3 elds six significant figurejs with a 

i 

maximum enor of three units m the seventh place Notice 

> 

that m the above discussion, some of the zeros in the second 

i 

and third components are counted as significant figures 
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One can similarly compute approximate values of y_ 2 (100.) 
and y 3 (100) 


(51) 


2 2 0O 



n=o 



"0“ 


' 0“ 


"-2 

po 

0 

1! 

1 

0 

* *1 = 

-1 

0 

’ —2 = 

-3 5 
-1 


a 3 


2 

3 8333*-- 
_1 


(52) 


Z 1 C100) 


532247 

X 

10 41 

2660306 

X 

10 45 

266124 

X 

10 41 


Comparing this value with the exact value, one can determine 

the accuracy to be seven significant figures with a three- 

unit error m the eighth place The actual error is roughly 
3 7 

3 x 10 which is a very large number This illustrates 
the situation where two large numbers can have several 
significant figures of agreement and still differ by a very 
large number. 

Similarly 


(S3) 


z 3 o) 



1 
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r 

— o 



" Cf 


3" 

0 

’ r 0 " 

-9 

-10 

— z 

49 


r 


3 


30 

90 

„-156 2/ 3_ 


(54) 


X 3 (100) 


- 100442 x 10~* D 

- 301326 x 10“ 45 
3365713 x 10 -42 


The exact solution yields the following information The 
approximate value (54) is accurate to five significant figures 
with an error less than two units m the sixth place 


For most problems there is of course no known exact (closed- 
form) solution Thus, one must attempt to determine the 
possible error involved by a more detailed investigation of 
the asymptotic expansions One technique is to merely 
truncate the expansion at the smallest teim (see [6]) Olver 
has recently investigated the problem of determining error 
estimates [9] anc it appears that, for certain cases, this 
procedure is quite bad, l e , the errors are much larger 
than one would expect This behaviour is most pronounced 
when one attempts to use the asymptotic expansions near a 
Stokes' line For the problem considered m this paper, 
the Stokes' lines correspond to the rays arg x = ± ^ 

Since the relative dominance of one solution over another 
changes when these lines are crossed, one might expect the 
inaccuracies that arise when one attempts to use the 
truncated foimal expansions 
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Note that the discussion so fdr depends on the relative 
dominance of one solution over another In fact, it is 
clearly possible to add any multiple of y (x) and £ 3 (x) 
to y_ 2 (x) without changing the asymptotic properties of 
y 2 (x) One could also add any multiple of y 3 (x) to 
£^(x) without changing the asymptotic properties of y^(x) 
Of course one must impose the lestriction | arg xj < ^ 

If one is interested m determining £ (x) and one only 
demands 


(55) 


CO 

y. (x) ~ / p x~ (n+1) where p 
n=o 



then the possible £ (x) form a one-parameter family of 
solutions {£ } Similarly, £ 2 (x) form a two-parameter 
family {£ 2 > Since no multiple of £ (x) and £ 2 00 can 
be added to y^ (x) and leave its asymptotic properties 
unchanged, one can conclude that this solution is uniquely 


defined 

(r - 

"0~ 

0 


r° 

1 


i.e , there corresponds precisely 


one actual solution of (23) asymptotic to the formal solu- 
tion 


NUMERICAL INVESTIGATION OF THOME' SOLUTIONS 

The following is a numerical procedure wh’ch, to* the author's 
knowledge, has never been published 

The only accurate information obtained concerning the solu- 
tions £ x ( x ) > £ 2 (x) and y 3 (x) is for large values of 
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X 


For many practical problems one is interested m these 

solutions for fairly small values of x In principle, one 

can use the values £.j£a) > £ 2 ( a ) and £ 3 (a) to continue 

these solutions back to small values of x Because of the 

different exponential rates of growth, tha s corresponds to 

✓ 

an unstable numerical problem Indeed, if one starts at 
x = 100 and attempts to use £^(100) to determine ££20) 3 
one will find that £ 3 (x) grows so rapidly (for numerical 
integration with a negative step the £ (x) solution of (23) 
dominates) that one soon obtains only a multiple of £ 3 (x) 

In ordex to obtain a vector £ (x) at x = 20 which is a 
member of the one-parameter family {£ } > one roust either 
increase the precision of the numerical calculations or 
develop some other numerical technique 

From the asymptotic relations one might expect 

(56) | |£ (x) | | ~ 05 at x = 20 

The growth of £ (x) is expected to be ex 10 34 where e 
is a sum of the errors m the initial condition (£^(100)) , 

the roundoff error, and truncation error Hence, one would 
expect to require about 40 places of accuracy to obtain four 
significant figures m the determination of £(20) . One 
would also require an extremely small step size m the 
calculations to maintain this precision 

Fortunately one can avoid this high-precision arithmetic 
One need only exploit the property that £ (x) can be 
considered a member of the one -parameter family {y } 
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Recall that any multiple of £ 3 (x) added to £ (x) yields 
a. solunon of (23) or (27) which is asymptotic to__ y^tx) > 
i e , the asymptotic properties are unchanged If one pro- 
ceeds to integrate from "x = 100 with a negative step, one 
will soon observe very rapid growth of the third component 
of the vector solution 


The rapid growth of the third component occurs since - £ (x) 
grows most rapidly m the third component If one CORRECTS 
£ (x) by adding a multiple of £ 3 (x) so as to force the 
third component to be zero at say x = 99 , then one has 
not changed the asymptotic properties of y 1 (x) One is 
merely calculating a different member (after the correction) 
of. the one-parameter family {£ } . If one repeats this 
procedure at x = 98 , x = 97 , • * • , x = 20 one arrives 
at a solution of the form 


(57) 


£^ 20 ) 


a_ 

i 

L° . 


Since there is also an exact solution for this particular 
problem one can determine the accuracy of this procedure 
(See Tables I and III ) 
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(NT) — numerically integrated 

33 

































For more practxcal problems there is no known closed-form 
solution Indeed, there would be no need to integrate 
numerically m such a problem One can attempt to calculate 
a solution asymptotic to y 2 (x) at x = 20 by numerically 
integrating from x = 100 to x - 20 For this problem, 
one adds multiples of (x) and X 3 ( x ) to correct y 2 (x) 

at x = 99 , x = 98 , •♦* , x = 20 At each of these 

values of x , one corrects y (x) so that it is of the 
form 


(58) 


l 2 (X) * 


‘0 

M x > 

0 


The value of g 2 so obtained at x = 20 can then be used 
to numerically integrate (23) with a positive step One can 
expect to maintain good accuracy on y (x) for integration 
with a positive step since £_ (x) dominates the other solu- 
tions of (23) m this direction When one reaches the value 
x = 100 , one can compare the value of (x) calculated 


using 


0 

3 2 ( 20 ) 
_0 


as an initial condition at x = 20 with 


the value one started the calculation at x = 100 These 
two values will not compare well unless one corrects the 
values of both to the same form (58) There is no reason to 


expect the initial condition 


0 

3 2 (20) 

0 


to he of the form 


(58) at x = 100 since there were 80 separate corrections 
and hence 80 different members of the two-pai ameter family 
{£ 2 > which have been used m the calculations The agree- 
ment (oi lack of it) at x = 100 is some indication of the 
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accuracy of all the numerical solutions obtained at x = 20 
This follows since all the solutions are used m the calcula- 
tion of y_ 2 (20) One may find that more precision is 
required, but it appears that one will require less precision 
than a numerical integration without utilizing corrections 
(see Tables II, III and IV) 

All numerical integrations were carried out using the UNIVAC 
1108, single piecision (eight significant figures) with a 
fifth-order Runge-Kutta method [7] and a fixed step size H 
(|H| = 1/32) 
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TABLE II 


X 

Z 3 (EXACT) 

z 3 (Mir 

100 


00010090 * 
- 00030278 
. 33657280. 

x 10" 42 


-- 0001004“ 
- 0003013 
3365713. 

x 10' 22 

90 


•- 00025738“ 
- 00077239 
69551972. 

x 10~ 38 


00025738' 
- 00077238 
69551542. 

x 10' 18 

80 


00006670“ 
- 00020017 
14244000. 

x ID' 33 


00006670“ 
- 00020017 
14243887. 

x 10" 13 

70 


*- 00017624' 
- 00052901 
28827363. 

x 10~ 29 


00017524' 
- 00052901 
28827086. 

x 10~ 9 

60 


00047740“ 
- 00143325 
57398994. 

x 10~ 25 


- 00047739“ 

- 00143323 
57398345. 

x 10" 5 

50 


00013357' 
- 00040113 
11161680. 

x 10” 20 


00013356' 
- 00040112 
11161535. 


40 


00039040“ 
- 00117316 
20912637. 

x 10' 16 


- 00039040“ 

- 00117314 
20912327. 

x 10 4 

30 


00121337“ 
- 00365093 
36685010. 

X lo” 12 


“- 00121335“ 
- 00365087 
36684401. 

x 10 8 

20 


- 00411079" 

- 0 1 241501 
53775284. 

X 1Q“ 8 


“- 00411072“ 
- 01241479 
55774266. 

x 10 12 


* 

(NI) — numerically integrated 


Notice the numerically integrated (NI) solution has been 
scaled No corrections are required to maintain good 
accuiacy on y (x) The dominance of y (x) for Integra- 
tion with a negative step completely explains tms accuracy 
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TABLE III 


X 

t x (EXACT) 

z x fNivcr 

^CNIW/OC)* 5 

100 


98010068' 
- 00009771 
_ 00029399_ 

X 10 2 

1 

BB 

X lo' 2 


98 010070' 
- 00009800 
.- 00029400 . 

x 10" 2 

90 

1 


X 10' 1 



X 10' 1 



x 10" 1 

80 


' 12189470' 
- 00001903 
00000000. 

X 10 _1 


12189426* 
- 00001903 
. 00000000. 

X 10' 1 


9 

x 10' 1 

70 


' 1388 04 95* 
- 00002830 
. 00000000. 

x 10 _1 



X 10” 1 


9 

x 10' 1 

BO 



x 10" 1 



X 10' 1 


16235352' 
- 00488281 
. OOOOOOOO. 

x 10' 1 

SO 



x 10 -1 



X JO -1 




40 

1 


x IQ" 1 


I 

X 10" 1 


9 

x 10 6 

30 



x 10 _1 



x 10 -1 


" 42949673" 
85899346 
. 00000000. 

x lO 10 

20 



x 10' 1 



X 10 _1 



x 10 16 


a (NIWC) — numerically integrated Mth corrections 
b (NIVI/OC) — numerically integrated without corrections 
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Notice for x = 90 , x = 80 , • , x = 20 all solutions 

are reduced to the same form 

) 

a (he)" 

P-lOO 
0 _ 


(59) 


ZiCx) - 


In order to facilitate a comparison one must examine the 
same member of the one -par ameter family {y } Thus, all 
vectors are reduced to the form (59) For the numerical 
integration with corrections (NIWC) , one uses this corrected 
value m the numerical integration, but the corrected value 
is ignored in the numerical integration without corrections 
(NIW/0C) Table III clearly indicates the value of the 
correction procedure 

One may wonder why the comparison of y (EXACT) and V (NI) 
m Table I seems to be much worse than y (EXACT) and 
y (NIIV/OC) The reason foi this is simply that one is 
obtaining a member of {y.^} when one computes y^(NI) 

One would need much more precision to more accurately 
specify a particular member of this familv at x = 100 m 

order to achieve agreement between y^ (EXACT) and y~(NI) 
This emphasizes the numerical instability of integration 
without corrections. 

* 

It appears that numerical integration without correction 
(NIW/OC) yields a member of the one-parameter family {y } 
to two significant figures at x = 60 This exceeds our 
expectations The reason is rather simple It jis expected 
that the first component of y (x) should be of the order 

_ 2 ^ i 

10 Since the inaccuiacies m the numerical calculations 

are going to be a result of rapid growth m the third com- 

< 
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ponent, one may determine how laige the third component will 
have to be m order to affect the first significant figure 
of £ 3 (x) From Table II one observes that £ 3 (be) is 
approximately of the foim 


(60) 




~-10 a 
-10 a 
- 1 0 a+3 _ 


Now, one expects to lose all significant figures m the 
determination of £ (x) when a = 6 (112^0011 ~ ) 

since one maintains only eight significant figures m the 
calculations Since the error m the initial condition is 
of the order 10~ 9 and the third component will have to 

9 

attain the approximate value 10 , this implies that a 

growth of 10 times the initial error is necessary to 

lose all significant figures The £ (x) solution of (23) 

18 ^ 

attains a growth of 10 times its original value for a 
decrease of x between 40 and 50 units Thus, at x = 60 
y^CNIW/OC) is a two significant figure member of {£ } 
and at x = 50 has no significant figures of agreement 
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TABLE IV 


X 

t 2 (EXACT) 

^(NlhC) 3 

| 

Z 2 CNI) c 

100 


‘ 00005319’ 
26603058 
. 00002658. 

< 10 45 

1 


X 10 25 

1 


x 10 25 


9 ^ 

x 10 25 

90 


0000000 0' 
11444181 
00000000. 

x 10 41 

1 


X 10 21 


■ 00000000’ 
11444162 
00000000. 

x 10 21 


• 00017569’ 
01208596 
S7 7 9224 6. 

X 10 22 

80 


' 00000000' 
48 9106S9 

. oooooooa 

x 10 36 

1 


x 10 16 


■ 00000000 
11197289 

oooooooa 

x 10 17 


‘ 00005542” 
00016633 
11835600. 

x 10 27 

70 



x 10 32 

1 


X 10 12 


■ OOOOOOOO” 
18447855 
oooooooo. 

x 10 21 


00014644" 

000439S7 

259»3141. 

X 10 31 

60 


’ 00000000’ 
86900848 
. 00000000. 

x 10 27 


■ 00000000’ 
86899977 
00000000. 

x 10 7 


■ oooooooo” 
90896679 
. OOOOOOOO. 

X 10 22 


' 00059668’ 
00119091 
47693709. 

x 10 35 

50 


' 00000000’ 
35878055 
. 00000000. 

x 10 23 

1 

jWy 

x 10 3 


NUMBERS TOO 
LARGE EXCEED 
10 38 

40 



x 10 19 

1 

rlSg it £ k HM 

x 10 _1 

|^S| 

x 

30 


" 00000000“ 
56362816 
. OOOOOOOQ_ 

x 10 14 

1 


X 10" 6 

X 


20 


■ 0000 0000" 
20434713 
. 00000000. 

x 10 10 



x 1 Q‘ 10 

ZxZ 

x 


®CNIWC) — numerically integral d with corrections 
b (NIW/OC} — numerically integrated without corrections 
e (Nl) — numerically integrated 
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From Table IV one concludes that £ 2 (NIW/OC) is much more 
accurate than y 2 (N I ) Since £ £NI) is used to compute 

y 2 (NIW/OC) one may wonder why there is such a discrepancy 
Before performing the comparison, one must put the solutions 
m the corrected foxm (58) If this is not done, one will 
be comparing different members of the two -parameter family 

<* a > 

One can predict the inaccuracies of £ 2 (NIW/OC) Since 
£ 2 (x) is decreasing exponentially fast for integration with 
a negative step, one expects a breakdown m the numerical 
solution, i e , no significant figure agreement with a 
member of {£ } for x near the value x = 80 

It is easily seen from the tables that for numerical integra- 
tion with corrections one obtains £ (x) to five significant 
figures with a seven unit error m the sixth significant 
figure, y 2 ( x ) t0 five significant figures with a four unit 
error m the sixth significant figure and £ (x) to four 
significant figures with a one unit error m the fifth 
significant figure 

If one uses the value of £ (NIWC) at x = 20 as an initial 
condition and integrates with a positive step to x - 100 

a 2 (100) 

then one obtains P 2 (100) If this value is 'reduced to 

Y 2 (100) j 

the form (58) and if the initial vector is also reduced to 

1 

this form one can compare these values The results of 
such a comparison yield four significant figure agieement 
with a one unit difference m the fifth significant figure 
This difference 3S seen to be of the same order as the error 

t 


41 



in y 2 (x) at x " 20 Since ^(x) and 2. 3 C X 3 aie used 
m determining this value of y_ 2 (x) , one might expect this 
difference to also be an indication of the errors present m 
these solutions 

It is lather surprising to obtain a better value of y_ 2 (x) 
at x = 20 than was obtained for y^(x) Apparently this 
is due to the relative inaccuracy of ^(100) as compared 
to y_(100) The inaccuracies m y (x) and y (x) 
slightly change the multiples of these solutions required 
for the correction process but do not seem to affect greatly 
the accuracy of y_ 2 (x) This is a very useful property 
since for larger systems, i e , for higher order systems 
one can hope to handle all the solutions m the above 
manner If the accuracy decreased, one would require 
increased piecision for higher order systems 

OBSERVATIONS AND FURTHER RESEARCH 

It is interesting to note that if one introduces the inner 
product 

C 61 ) 

where y_ 
the angle 

(62) 


Cxu) ■ + y 2 v 2 - y 3 V : 


1 

X 

H 

1 


" V l" 

y 2 

and v - 

V 2 

— i 

ro 

>-. 

- i 


V 3 


then one can introduce 


9 , 


0 = arc cos 


(X>v) 


/Ty yZ) ” x ~ (X>X3 
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The corrected solutions £ (x) > £ 2 ^ > £ 3 (x) aie practi- 
cally mutually orthogonal The minimum angle occurring for 
any pair of solutions (100 <_ x < 20] is nearly 89° The 
correction process is roughly equivalent to a Gram-Schmidt 
process Since a method of this type has already been 
investigated by Conte [5] it would be interesting to see 
the method m this paper and Conte's method compared In 
this connection, it is even more interesting to note that 
this author has more nearly satisfied the conditions that 
Conte wanted to satisfy, l e , Conte was willing to tolerate 
1° values of 6 

One should^ exploit, if possible, any freedom that a problem 
allows For example, one can define an inner product so 
that an arbitrary set of n linearly independent vectors 
in an n- dimensional space form an orthonormal basis One 
should choose this basis after one has simplified the first 
few coefficient matrices [8] , 

Since Conte's method and the method m this paper m th.eir 
present forms seem to be intended for use in different 
problems one would first need to extend both methods before 
attempting a comparison 

If one attempts to utilize the correction procedure 
described m this paper one must be careful to integrate m 
a direction so that the correction procedure is justified 
For example if one attempted to integrate with a positive 
step and obtain £ (x) one wou ld anticipate the need to 
correct \\fith a multiple of £ 2 ( x ) Since £ 2 (x) dominates 
£ (x) asymptotically there is no justification for adding 
multiples of y (x) to £ (x) , i e , one does not mam- 
tain the same asymptotic propeities 
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